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Abstract 

In this paper, we consider a Levy-driven fluid queueing system where the server may 
subject to breakdowns and repairs. In addition, the server will leave for a vacation each 
time when he finds an empty system. We cast the queueing process as a Levy process 
modified to have random jumps at two classes of stopping times. By using the Kella- 
Whitt martingale method, we obtain the limiting distribution of the virtual waiting 
time process. Moreover, we investigate the busy period, the correlation structure 
and the stochastic decomposition properties. These results may be generalized to 
Levy processes with multi-class jump inputs or Levy-driven queues with multiple input 
classes. 

Keywords: Levy processes; Fluid queues; Server breakdowns and vacations; Kella-Whitt martin¬ 
gale; Stochastic decomposition 

Mathematics Subject Classification: 60K25; 60K30 

*E-mail address: wujinbiao@ymail.com 


1 



1. Introduction 


The class of Levy processes consisting of all stochastic processes with stationary independent 
increments is one of the most important family of stochastic processes arising in many areas of 
applied probability. It covers well-studied processes such as Brownian motion and compound 
Poisson processes. Levy processes are used as models in the study of queueing systems, insurance 
risks, storage systems, mathematical finance and so on. For a comprehensive and outstanding 
analysis of Levy processes and their applications, readers can refer to the books by Bertoin 14fl, 

Q. - n 


Sato 


22], Applebaum [if] , and Kyprianou 1181- 


Over recent years it has been a rapid growth in the literature on queues with server breakdowns 
and vacations due to their widely applications in computer communication networks and manu¬ 
facturing systems. Li et al. [19] considered an M/G/l queue with Bernoulli vacations and server 
breakdowns using a supplementary variable method. Gray et al. Q analyzed a multiple-vacation 
M/M/1 queueing model, where the service station is subject to breakdown while in operation. Ke 
151 ] studied the control policy in M^ /M/1 queue with server breakdowns and multiple vacations. 


iin [14] 


Jain and Jain M dealt with a single server working vacation queueing model with multiple types 
of server breakdowns. Wu and Yin Q gave a detailed analysis on an M/G/l retrial queue with 


non-exhaustive random vacations and unreliable server. Wang et al. 


231 ] studied a discrete-time 


Geo/G/1 queue in which the server operates multiple vacations and may break down while working. 
Yang and Wu [25] investigated the N-policy M/M/1 queueing system with working vacation and 
server breakdowns. However, most of the papers are focused on queues where customers arrive at 
the system according to independent Poisson processes or geometrical arrival processes. Queue¬ 
ing systems with Levy input (or Levy-driven queues) are still not well investigated in literature. 
Recently, queuing systems with Levy input have been attracting increasing attention in the ap¬ 
plied probability and stochastic operations research communities. Levy-driven queues covers the 


classical M/G/l queue and the reflected Brownian motion as special cases. Kella and Boxma 


Q 


first considered Levy processes with secondary jump input which were applied to analyze queues 


with server vacations. Recently, Lieshout and Mandjes 


20] analyzed tail asymptotics of a two-node 


tandem queue with spectrally positive Levy input. Boxma et al. [5] analyzed a generic class of 


• S 


Levy-driven queuing systems with server vacation. D§bicki and Mandjes [9| provided a survey 


on Levy-driven queues. D§bicki et al. [8] focused on transient analysis of Levy-driven tandem 


2 










queues. Palmowski et al. 


fl 


21[ considered a controlled fluid queuing model with Levy input. Boxrna 


and Kella [6|] generalized known workload decomposition results for Levy queues with secondary 
jump inputs and queues with server vacations or service interruptions. However, to the best of our 
knowledge, no work on Levy queues with server breakdowns and vacations is found in the queueing 
literature. 


In this paper, we consider a single-server Levy-driven fluid queue with multiple vacations (ex¬ 
haustive service) and a server subject to breakdowns and repairs which is motivated by the perfor¬ 
mance analysis of resources in communication networks. A breakdown at the server is represented 
by the arrival of a failure. The principal purpose of the present paper is to apply the martingale 
results which were derived in Kella and Whitt [f/J to investigate the stochastic dynamics of the 
system and realize an extensive analysis of the system from the transient virtual waiting time pro¬ 
cess to the steady-state distribution of the waiting time process. In addition, we give a detailed 
analysis of the system busy period and the queue’s correlation structure. Furthermore, we establish 
decomposition results for the model. 

The rest of the paper is organized as follows. In Section 2, we introduce a few basic facts 
concerning Levy processes and martingales. In Section 3, we describe the general model and 
formulate the model as a Levy process modified to have random jumps at two classes of stopping 
times. In Section 4, we characterize the steady-state distribution of the virtual waiting time process 
and the mean length of the busy period. In Section 5, we address the transient distribution of the 
waiting time process as well as the queue’s correlation structure. In Section 6, we present two 
stochastic decomposition results. 


2. Preliminaries on Levy processes 


In this section, we mainlygive several notations, definitions and propositions about Levy pro- 


rr. 

cesses. For more details, see [1|, 2, 4, 


22 ], 


Denote M = (— 00 , 00 ) and M+ = [0, 00 ). Let (D, F, P, {^> 0 }) be a complete probability space 
endowed with a cr-field filtration Ft> 0 , he. an increasing family of sub-fields, which fulfils the usual 
conditions. That is, each Ft is P-complete and Ft = n s> tF s for every t. Throughout, adapted, 
stopping times and martingales will be defined with respect to this filtration. 
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Definition 2.1. Let X = (X t ,t > 0) be a stochastic process defined on a probability space (fl, J 7 , P), 

if 


(i) X$ = 0 (a.s.); 

(ii) X has independent and stationary increments; 

(Hi) X is stochastically continuous, i.e. for all a > 0 and for all s > 0 ; 

limP(|X t — X s \ > a) = 0, 

t—ys 

then X is a Levy process. 

Two central and key results in the foundations of Levy processes are given by the following 
propositions (known as the Levy-Khintchine formula and Levy-Ito decomposition respectively). 

Proposition 2.1. Let X be a Levy process with Levy measure v (J K _| 0 j.(|x| 2 A l)v(dx) < oo). 
Then 

Ke idXt = e^ (0) , t> 0, 

where 

ip(6) = ifiO — + f (e ldx — 1 — i9xl^ < i)n(dx), 

2 J R-{0} 

/j€M and a > 0. Moreover, given u, /i, a 2 , the corresponding Levy process is unique in distribution. 
Furthermore, the jump process of X is a Poisson point process with characteristic measure v, that 
is 

n(A) = E ]T 1 a (AX s ) . 

0<s<l 

The function 'if(-) is called the Levy exponent of the Levy process. 

Proposition 2.2. Let X be a Levy process. Then X has a decomposition 

X t = B t + p.t+ / x(N t {-,dx) - tu(dx)) + / xN t (-,dx ) 

J{\x\<l} J{\x\>i} 

where B is a Brownian motion and for 0 ^ A, N = f A N t (•, dx) is a Poisson process independent 
of B. 

The following are some known fundamental results about Levy processes: 
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(i) Any Levy process can be represented as an independent sum of a Brownian motion B and a 
‘compound Poisson’-like process A" 0 . If A has no negative jumps and the paths of A 0 are of 
bounded variation, then A 0 can be a subordinator which can be represented as a nonnegative 


compound 
processes 


'oisson process or as the limit of a sequence of nonnegative compound Poisson 


lall- 


(ii) If z/((—oo,0)) = 0, then the Laplace-Stieltjes transform exists and is given by Eexp(— OXt) = 
exp (t<p(9)) where 

<p(0) = logEexp{ — 9X{\ = — /j,0 + ^-0 2 a 2 + f (e~ dx — 1 + 0xl\ x \ < i)v{dx), 

2 J R + 

here the function </?(■) is called the Laplace-Stieltjes exponent of the Levy process. It is easy 
get that 99 ( 0 ) = 0 and <£>(•) is convex; If X is not a subordinator, then ip(0) -> 00 as 0 -> 00 ; 
Whenever EX\ < 0, ip(-) is strictly increasing. Throughout the paper, exponent will mean 
the Laplace-Stieltjes exponent ip(-) due to the convenience compared to the Levy exponent 
</>(•)■ 


(iii) If A has bounded jumps, then E|Ai| n < 00 for every n > 0. 


It is known that there always exists a version with sample paths of A that are cadlag, which is 
therefore strong Markov. Throughout this paper, every Levy process mentioned is assumed to be 
such. 


In this paper, we focus only on the spectrally positive Levy process, i.e. v{— 00 ,0) = 0. For 
any spectrally positive Levy process A with (3 > 0 and EX\ < 0. Denote by T^ = inf-{> | At = -£} 
for any non-negative random variable £ which is independent of A. Then it follows Theorem 3.12 
in [ij| that 

£ex p{-/?T«} = expl-^CS)?}, ET^ = -|^-. 


We now give the generalized Pollaczek-Khinchine formula associated with a reflected Levy 
process with no negative jumps which is also a celebrated formula in queueing theory. Given some 
random variable 0 < Rq S Eq, let 

It = (— inf X s — -Ro^) , Rt = Ro + Xt + It, t > 0 . 

V 0 <s<t J 

Then R is a reflected Levy process with I being its local time at zero. Since A has no negative 
jumps, I is continuous with Iq = 0 and is the minimal right-continuous nondecreasing process such 
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that Rt > 0 for all t. The process {It, t > 0} must only increase when R t = 0, so that 


l{i? t >o}^ t — 0. 


The following result is the famous (generalized) Pollaczek-Khinchine formula in queueing theory. 
We can refer to the papers {3] and [ii| for its proof. 


Proposition 2.3. If X is a spectrally positive Levy process such that EX\ < 0, then 


lim Ee~ 9Rt 

t—t oo 


V( 0) 

m 1 


e > o. 


We end this section by presenting the famous Kella-Whitt martingale associated with spectrally 
positive Levy process X with exponent p(9) that we will apply it to analyze our fluid queue. 

Proposition 2.4. Let X be a spectrally positive Levy process with exponent ip(9). Let L t = J'q dL^ + 
Ylo< s <t A L s be an adapted cadlag process of bounded variation on finite intervals with continuous 
part {Lf} and jumps A L s = L s — L s _, and define Z t = Zq + X t + L t , where Zq € Fq. Then 

exp {—9Z s }ds + exp{— 9Zq} — exp {—9Z t } — 9 [ exp {—9Z s }dYf 

Jo 

+ ^2 [exp {-9Z S } - exp {-9Z S - AL S }] 

0 <s<t 

is a local martingale. In addition, if the expected variation of {L},t > 0} and the expected number 
of jumps of {L t ,t > 0} are finite on every finite interval and {Z t ,t > 0} is a non-negative process 
(or bounded below), then {M t ,t > 0 } is a martingale. 



The Kella-Whitt martingale is a stochastic integral with respect to the Wald martingale Nt = 
exp{— QXt — ip(9)t} and the readers may refer to the paper Q for its proof. 


3. Model formulation 

We consider a single fluid queue with a (nondecreasing) Levy input and available-processing 
(or service) processes. For any t > 0, let Xt be the cumulative input of fluid over the interval [0. t]. 
We assume that A is a subordinator, that is, a Levy process, having the exponent 

<f>{9) = — a9 + f ( e~ dx — 1 )v{dx), 

J R+ 
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where a > 0 and f R xv(dx ) < oo. Let S(t) be the cumulative available processing over the interval 
[0, t\. Here, we assume the server processes the fluid at a constant deterministic rate r (whenever 
the fluid level is positive). So that S(t) = rt, t > 0. We assume the storage space is unlimited. 
Let {r n , n > 1} be the strictly increasing sequence of stopping times at which the server fails. 
When the server fails, it is repaired immediately and the time required to repair it is a positive 
random variable which is -measurable for n > 1. When the system becomes empty, the 
server takes a vacation of random length r/ n which is J - an -measurable for n > 1, where {cr n ,n > 1} 
is a strictly increasing sequence of stopping times at which the server is leaving for a vacation. 
Vacations continue until, on return from a vacation, the server finds the system is non-empty. 



Remark 3.1. It should be noted that here Wt, t > 0, is not the workload process. However, if 
represents the total arriving workload during the nth repair period and p n corresponds to the total 
arriving workload during the nth vacation period, then Wt is the workload in the system conditioned 
on being in the active periods which are defined as the normal processing fluid periods (excluding 
the repair periods and vacation periods). 


Define Y) = Xf—rt. Then Y is a spectrally positive Levy process with exponent y>{0) = (j>(0)+rO. 
Denote p = — ^(0) and impose the conditions p < r and E|Y)| < oo for all t > 0. 


4. The steady-state distribution 

In this section, we characterize the limiting distribution of the virtual waiting time Wt- Let 
—>• denote convergence in distribution and let —* denote convergence in probability. Under the 
condition p < r, we suppose Wt A W and t -1 E Wt —>• 0, as t —>• oo, where W is a random variable. 
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In order to derive the steady-state distribution of the virtual waiting time Wt, we give the 
following lemmas. 


Lemma 4.1. If E(N^ + ) < oo for all t, then {Mt,t > 0} is a zero-mean real-valued martingale 

with respect to {Ft, t > 0 }, where 


M t = tp(0) [ e~ eWs ds + e~ ew ° - e~ 9Wt 

Jo 


J2 - e~ 0W ^ 

k= 1 


Nf 

E( i - *-**)■ ( 4 -‘) 


isjR mV 

Proof. Let H t = Y2i=i + XAi Vi- Then Wt = Wo + Y t + H t and {H t , t >} is adapted of bounded 
variation on every finite interval, and H£ = 0 for all t > 0. Since E{N^-\-NY) < oo, by Proposition 
12.41 and through a trivial manipulation, we get that H4.ip is a martingale. □ 

Lemma 4.2. Let T is a nonnegative random variable with possible integer values 0,1,2, ••• and 
E T < oo. If {C n ,n > 1} is an i.i.d. sequence with independent of 1 {r>n} an< J ECi < oo, then 

T 

Ej]Cn = (ECi)(ET). 

n =1 

Proof. Since f n is independent of 1 {T>n}, the bounded convergence theorem yields 

T oo oo oo 

Cn] = E[y^ Cnl{T>n}] = ^(^Cn)(!El{T>n}) = ®Cl y~](El{ r > n |) 

Tl= 1 72—1 72 — 1 72 — 1 

OO 

= ECi E P(T >n) = (ECi)(ET). 

72—1 

□ 

With the help of the above two lemmas, we can derive the limiting distribution of the virtual 
waiting time provided in the following theorem which is the most important result of this paper. 


Theorem 4.1. Let {£ n ,n > 1} and {? 7 n ,n > 1} be two positive i.i.d. sequences with E£i < oo and 
Er/i < oo, respectively. Suppose that n~ l r n A A ^ 1 and n~ l a n A Ay 1 as n —» oo for 0 < A ^ 1 < oo 
and 0 < Ay 1 < oo, respectively. Further, assume that {W Tk , k > 1} and {W Tk — £fcA > 1} are 
stationary and ergodic that (W Tk — £k,W Tk ) A [W~ ,W + ) where W~ ,W + are two proper random 
variables. If is independent ofl^ N R >n y and r] n is independent of l^ N v >n y, then for 6 > 0, 


_ pA( o) x (i -p)A(Q) 

E£i ’ ' Et7! 


(4.2) 
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(4.3) 


N R 


N* 


lim t £*. = —pEY\, lim t rik = — (1 — p)Eli, 

t —loo ^ £—loo A—^ 


fc=l 


£—>•00 


fc=l 


lim Ee- 0H/ * = Ee- 0H/ = 

<—>■00 ^(p) 


' Ee- 0W " - Ee- 0W+ 

P — 


0 E& 


N 1 - Ee _07?1 
+ (1 ~ i,) 0 E, a 


where 0 < p < 1. 


(4.4) 


Proof. Dividing (13.11) by t and taking expectation, we get 


N? 


nY 


t^EWt = t~ l EW\ + Eli + t _1 E ^ %. (4.5) 

fc=i fe=i 

where t _1 EHi —>• 0 by the assumption. Thus, we have proved ([4.31b 

Taking Wq = W* where W* is an independent copy of W, {Wt,t > 0} becomes stationary. By 
Lemma 14.21 and (13.11) , we get 


ElVfE^i + E Nf E 771 = -Eli. 


Vn 


(4.6) 


So that, 


ElVfE^i = -pEYi , ENf Er/! = - (1 - p)EY x . 


Vv 


Hence, 


R _ PP'i °) w AT V _ U “ pW (°) 


ElVf = 




EiVj = 


Er/i 


(4.7) 


Conditions n 1 r n A A ^ 1 and n 1 a n A Ay 1 as n —>• 00 imply that t 1 Nf' Aa« and t 1 iN/^ 
Ay. Together with t~ 1 Nf t —> EA^ and t~ x N^ —> EA r ) w.p.l as £ —>• 00 , this implies (14.21) . 

From Lemma 14.21 optional stopping at t = 1 in (14.11) yields 


0 = ^(6»)Ee" eM/ - E1V/ 


ri? 


Ee 


- 6 W~ 


— Ee 


-ew+ 


-EA^l-Ee" 0?71 ) 


Substituting ([4.7D into (14.81) . we obtain (14.41) . 


(4.8) 


□ 


Remark 4.1. From \\], we have <p(9) = 0 for some 0 > 0 if and only ifY t has a lattice distribution. 
But for a Levy process with no negative jumps and Eli < 0, it is not possible. In addition, when 


9 












A T t R = 0 for all t, i.e., p = 0, the present system reduces to the Levy process-driven queue with server 
vacations. It may be noted that the equation {HD after putting p = 0 agrees with the equation (f.6) 
presented in Kella and Whitt [16], We should note that when = 0 for all t, the present system 
reduces to the Levy process-driven queue with server breakdowns. However, the equation (14.4|) in 
Theorem \f.l\ after putting p = 1 cannot characterize the limiting distribution of Levy-driven queue 
with server breakdowns due to the fact that {Wt,t > 0} is not a non-negative process. 


Corollary 4.1. When the system is stable, we have 


u = EW = f^ +p nW^-E(W-r + *4 


Mo) 


2E£i 


2IE?7i ' 


v =VarfT = 


1 ^"(0) 

3 7(0) 


^(0) 

^(o) 


+ (1 -p) 


1 Er/ 3 
3 E i/ i 


MV 

E771/ 


+ p 


1 E(1C+) 3 - E(iy-) 3 

3 E6 


1 /E(IT+) 2 -E(W~) 2 
4 V E^ 


Next, we analyze the system busy period which is defined as the time period which starts at the 
epoch when a vacation is completed and the server begins processing the fluid and ends at the next 
epoch when the system is empty. Note that a busy period includes the normal processing time of 
the fluid in the system and some possible repair times of the server due to the server breakdowns. 


Lemma 4.3. Define T n = inf{t|Yt + ^) fc= i 1 + -B* = 0}, n > 1, where {B- L ,i > 1} is a positive 

i.i.d. sequence with EBi < oo. If 0 < p < 1, then 


E T n 


nEB\ 

(i-pV(o)’ 


(4.9) 


Proof. Let = inf{t|F t + Ylk =i = 0}, b > 0, n > 1. Since Y is a spectrally positive 

Levy process, { T b , n > 1} is a random walk with 0 < T b < T b +1 w.p.l. Thus, E T b = nETfk Since 

,__ ]\jR 

E Y t < 0 and 0 < p < 1, E(Y* + Ylk =l £fc) = — (1 — p)p'{0)t. < 0. By Theorem 8.4.4 of Chung (1974), 

we have ET^ < oo. Since E|Y)| < oo, t~ l Yt —>• Eli w.p.l as t —>• oo by the Kolmogorov strong law 

N* b 

of large numbers. Taking into account the fact that Y T b + Ylk=i = ~ we get 



( \ 


1- 

fe; 

1_ 

b = n 1 

Y n + 

= n~ l T b n 



\ k=1 


k =1 


Erf[EYi + (EATf )(E^i)] = ET^[-(1 - p)<p'( 0)] w.p.l as n -> oo. 
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So that ET^ = • Finally, by conditioning and unconditioning, we obtain (|4.9h . □ 

By the argument of Lemma 14.31 we obtain the following result. 

Theorem 4.2. Let T = 'mi{t\Wt = 0}, where Wq is distributed according to the stationary distri¬ 
bution. Then the mean length of the busy period is given by 

EW 

ET = 


(i -pV(o)' 

Remark 4.2. Note that for the M/G/l queue with server breakdowns and multiple vacations, if 
Bi is distributed as a service time, then T n is the n-order busy period which is defined as the time 
period which starts when a vacation is completed and there are n customers in the system and ends 
at the next departure epoch when the system is empty. 


5. The transient distribution 


In this section, we focuses on analyzing the transient distribution of Wt in terms of Laplace- 
Stieltjes transform, for some t > 0, conditional on Wq = x. 

Theorem 5.1. Let T be exponentially distributed with mean 1/7, independently ofY, r n and cr n , 
n > 1. For 9 > 0 and x > 0, 

E x e~ eWT 


ip( 9 ) - 7 


Ee~ 0w ~ - Ee- 0W+ 


1 - Ee~ dm 


X i ^Ee - ^ -1 ^)^ - - Ee~v~ 1 G)W+ + ^ P \ - 6 


_ P -[8-p Ht)]® 


Proof. By Lemma 14.11 we have 


roc rt 

<p{6) / / ^e-^e- eWa dsdt + e~ 6x -E x e 

Jo Jo 


0 = EMt = 


-0W T 


N* 


-E^ - e~ ew ^ 

k=l 


Nf 


*E(‘ 


— e 


SVk' 


k =1 


ip( 6 ) 


7 


E xe -ew T + e -ex _ E _ ^ e -ew T _ EN R 


Ee 


—8W~ 


— Ee 


-ow+ 


-EN^(1-Ee~ er11 ). 


So that 


E x e~ eWT = —l - {e Nfi 

ip( 9 ) - 7 l 


Ee 


-ew~ 


— Ee 


-w+ 


+ EN%(1 - Ee~ 0r>1 ) - e~ 9x j . 


(5.1) 
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Note that when EYi < 0, ip{0) is increasing on [0, oo). Therefore, the inverse of ip(9 ) is well dehned 
on [0,oo). Hence, the equation ip(9) = 7 has exactly one root 9 = </? - 1 ( 7 ) on [0,oo). Finally, using 
the fact that the root of the denominator should be a root of the numerator as well, we can obtain 

ED- □ 


Next, we derive explicitly the Laplace transform corresponding to the correlation of the virtual 
waiting time process: 

r(A = Cov(W 0l W t ) = E(W 0 W t ) - (EWp ) 2 = E (W 0 W t ) - cn 2 

’ VVarVFo • YarW t YarW 0 v 

Here, we assume the system is in steady-state at time 0. 


Theorem 5.2. For 9 > 0 and u, v as in Corollary \4-l 


c{t)e et dt =/- — 


1 <V(0) E (We~v 1 ^ w ) 


8 v8 2 

x 


+ 


v9 


(1 ~p)Et?i 


pE^i 


+ 


where 


E (We-v 1(0)M/ ) = 


<y9'(0) 9ip'(0)<p'{9) 


.<P( 8 ) ip 2 (9) 

p'{ 0) P 


Ee 


—6W~ 


P- 


— Ee 


-6W+ 


+ 


+ 


p{8) ECi 

l'(0) 1 ~P 
<p(0) Er /1 


E (W + e~ ew+ ) - E (W~e- ew ~) - 
1 - Ee~ dm 


+ (1 ~P) 


1 - Ee ~ e ^ 


9Ern 

Ee~ ew ~ - Ee~ ew+ 

9 


E{ m e- eri1 ) - 


8 


Proof. Let T be exponentially distributed with mean 1/9. From (j 5.1 h . we get 


f 9e~ et E x W t dt = - + e-* 1{e)x 

Jo 9 


_ (1 -p)Ern 

_ E e -v 3_1 ( e )W /+ ' i_Ee _ ^ _1 ^)m 




+ 


Straightforward calculus yields 

/* OO I pOO 

/ c{t)e~ et dt = - [E (W 0 W t ) - 0 J 2 ]e~ dt dt 

Jo v Jo 


1 


v Jo Jo 

Substituting (|5.2D into m , we obtain 

,2 


poo poo 

/ / xE x W t e~ et dF(W 0 < x)dt - 

Jo Jo 


c 0 
v9 


poo 

/ c(t)f 
Jo 


e et dt = --+ 

v9 


uj 


+ X. 

(5.2) 


(5.3) 
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[°° x_ 

Jo vO 


e 


_|_ e -¥> 1 (P)x 


(1 -p)Erj 1 
- E e -‘P~ 1 (0) w+ ' l - Ee-v 3-1 ^)^ 1 


pMi 


+ 


+ x| 6F{Wq < x) 


1 cV( 0) E^ITe”^ 1 W ff ) 




+ 


(1 -p)Er?i 


■u6» 2 ' 1 Ec-r 1 ^- -Ee-r 1 ^ 1 - Ee-* 5-1 ^ 

where E(I / Le _¥3 can be obtained by differentiating (14.411 . This completes the proof. □ 


6. Stochastic decompositions 


In this section, we identify some stochastic decomposition properties of the virtual waiting 
time for our fluid queueing system. Stochastic decomposition properties have been studied in 
many vacation models. The classical stochastic decomposition property shows that the steady- 
state system size at an arbitrary point can be represented as the sum of two independent random 
variable, one of which is the system size of the corresponding standard queueing system without 
server vacations and the other random variable depends on the meaning of vacations in specific cases 
(see Fuhrmann and Cooper [11] and Doshi [1(J). In addition, stochastic decomposition properties 
have also been held for Levy-driven queues with interruptions as well as cadla g p rocesses with 
certain secondary jump inputs (see Kella and Whitt jl6 ] and Ivanovs and Kella 


13]). It’s worth 


mentioning that, recently, Boxma and Kella p] generalize known workload decomposition results for 
Levy queues with secondary jump inputs and queues with server vacations or service interruptions. 
In particular, in the context of our fluid queueing system, we have the following decomposition 
results. 


Lemma 6.1. There are two random variables U and V such that 


Ee 


-ew~ 


— Ee 


-ew+ 


mi 


= Ee 


-0U 


( 6 . 1 ) 


1 - Ee~ 6m 

OErji 


= Ee 


-ev 


( 6 . 2 ) 


Proof. The left hand side of (|6.1D is the Laplace transformation of the function 

P(W+ > x) - E(W~ > x) 


/(*) = 


Mi 


which is a bona fide probability density function. 
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The left hand side of (16.21) is the Laplace transformation of the function 


g(x) = 


P(r/i > x ) 

Ei]i 


which is a bona fide probability density function (stationary residual life density of 771 ). □ 

Applying Lemma l6.ll we can get the following theorem to characterize the stochastic decom¬ 
position property. 


Theorem 6.1. Under the conditions of Theorem \4-l\ the distribution of W is the convolution of 
two distributions one of which is the distribution of R in Section 2 and the other is the mixture 
of the distribution of U with probability p and the stationary residual life distribution of r/i with 
probability 1 — p. 


Furthermore, under some conditions, we have another stochastic decomposition result. 


Lemma 6.2. Under the conditions of Theorem \f.l\ if positive random variables f n , n > 1 are 


independent of Wo, Y Tri and Y2i=i Vi, then we have 

<V(0) T 1 -Ee "^ 1 


Ee 


-ew 


<p (6 


0E£i 


Ee 


-ew~ 


+ (1 -p) 


1 - Ee~ 6m 

@E?7i 




Proof. Since f k is independent of W Tk - f k = W 0 + Y Tn + JT' =1 fi + J2i= 1 Vi, 

Ee- d (W n -W _ Ee ~ ew n Ee~ e ^-^\l - e~ e ^) 


BEt! 


0E& 

= Ec~ e{WT ' ~ gl) 1 ~ Ee ~ ga 
1 - Ee”^ 1 


= Ee 


—9W~ 


OMi 


Substituting (16.41) into (14.41) yields (16.31) . 


(6.3) 


(6.4) 


□ 


From Lemma 16.21 we obtain the following theorem which provides another stochastic decom¬ 
position. 


Theorem 6.2. If the assumptions of Lemma 6.2 are satisfied, then the distribution of U is the 
convolution of the distribution ofW~ and the stationary residual life distribution off 1 . 
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